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Abstract 

We derive a Toda-type recurrence relation, in both high and low temperature regimes, 
for the A - extended diagonal correlation functions C{N, N; A) of the two-dimensional Ising 
model, using an earlier connection between diagonal form factor expansions and tau-functions 
within Painleve VI (PVI) theory, originally discovered by Jimbo and Miwa. This greatly 
simplifies the calculation of the diagonal correlation functions, particularly their A-extended 
counterparts. 

We also conjecture a closed form expression for the simplest off-diagonal case C*(0, 1; A) 
where a connection to PVI theory is not known. Combined with the results for diagonal cor- 
relations these give all the initial conditions required for the A-extended version of quadratic 
difference equations for the correlation functions discovered by McCoy, Perk and Wu. The 
results obtained here should provide a further potential algorithmic improvement in the 
A-extended case, and facilitate other developments. 
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1 Introduction 



The 2D Ising model in zero magnetic field is arguably the most important solvable model in 
theoretical statistical mechanics. Its free energy [1] and the spontaneous magnetization [2, 3] 
were computed a long time ago by Onsager and Yang. 

Despite more than half a century of attempts by many outstanding scientists, a closed form 
expression for the susceptibility is still unknown. We have a vast amount of numerical and 
analytical information, mostly derived from careful analysis of exact series expansions. The best 
algorithm for such series expansions, due to Orrick, Nickel, Guttmann and Perk [4], is based 
on expressing the susceptibility in terms of correlation functions, and using quadratic difference 
equations, due to McCoy, Perk and Wu [5], [6], [7] to recursively calculate the required correlation 
functions. A related study that gives complementary analytic information is based on the multi- 
particle expansion [8] of the susceptibility derived by Wu, McCoy Tracy and Barouch in 1976. 
In either formulation, one requires the correlation functions or their A- extended counterparts 
(see (7-8) below). 

In this paper we present a different method for the calculation of diagonal pair correlation 
functions and form factor expansions that occur in the theory of the 2D Ising model in zero 
magnetic field. 

We shall restrict ourselves to the case of the symmetric 2D Ising model defined by the 
interaction energy 

E = - J ^(o-jjo-jj+i + aijai+ij), (1) 

where J > (ferromagnetic regime) and spins aij = ±1 are assigned to each site of a square 
lattice. 

The 2D Ising model has a second order phase transition at the critical temperature Tc defined 
by the condition 

sinh(2J//cBrc) = 1, (2) 

where fcs is Boltzmann's constant. 

A study of the pair correlation functions 

C{M,N) = {ao,o(TM,N) (3) 

started long ago [9] when they were expressed in terms of determinants [10, 11]. Subsequent 
studies of correlation functions [7, 8, 12-17] led to infinite form factor expansions with every 
form factor given by a multiple integral [8]. 

Note also that there are two special cases where a simpler determinant representation for 
correlation functions is known. These are the cases of row and diagonal correlations M = and 
M = N respectively, for any N. In both cases there exists a representation for pair correlation 
functions hy N x N Toeplitz determinants [11]. 

In 1980 Jimbo and Miwa [18] found yet another approach to the calculation of the diagonal 
correlations C (N, N) . They showed that C {N, N) satisfies the "sigma" form of the Painleve VI 
equation and produced a system of recurrence relations which allows one to calculate diagonal 
correlations C{N + 1,A'" + 1) in terms of C{N,N) and C*{N,N) which are the correlation 
functions obtained by the duality transformation interchanging high/low temperature regimes. 

Moreover, Jimbo and Miwa introduced the so-called A-extension (or isomonodromic defor- 
mation) of the correlation functions using a multi-particle expansions in terms of free fermion 
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operators. Then they related this to the matrix Schlesinger equations with integrabihty con- 
ditions given by the Painlevc VI equation.^ In such an approach the parameter A describes 
monodromy properties of solutions for Schlesinger equations and plays the role of boundary 
conditions for the Painleve VI equation. 

Similar A-extensions naturally appear in the Predholm determinant approach [8] and can 
be regarded as generating functions for the form factors (see (7-8) below). The case A = 1 
corresponds to the symmetric Ising model. 

Starting with the asymptotics of diagonal correlations in the high/low temperature regimes, 
the authors of [20] studied series solutions of the Painleve VI equation and showed the existence 
of the one-parameter extension of C{N,N) which exactly matches the form factor extension. 
They conjectured that the Jimbo-Miwa A-extension and the form factor A-extension are actually 
the same, with the parameter A in (7-8) playing the role of boundary condition for the Painleve 
VI equation. 

We note that such an identification is a result of experimental mathematics, confirmed by a 
large amount of computer calcTilation. In this paper we will rely on this, even though a rigorous 
mathematical proof of such an identification is still missing. 

A different approach to calculation of correlation functions was developed by McCoy and 
Wu [5] and by Perk [7] in 1980. They discovered a set of quadratic difference equations satisfied 
by the pair correlation functions of the 2D Ising model. These equations were developed further 
in [21]. 

In the symmetric case the quadratic difference equations allow one to calculate all correlation 
functions recursively, once we know the "initial conditions" : diagonal correlations C (N, N) and 
C*{N, N) together with the nearest-neighbor correlations C(0, 1) and C*(0, 1) [22]. The diagonal 
correlations can be calculated from the Jimbo-Miwa recurrence relations and when combined 
with well known expressions for the nearest-neighbor correlations, gives a very effective method 
for series calculations of the correlation functions and the susceptibility [4] . 

In this paper we present a different method of calculation of diagonal correlations in the 
A-extended case. It is based on the explicit relationship of such correlations to tau-functions in 
Painleve VI theory. We also conjecture a formula for the A-extension of the nearest neighbor 
correlations C(0, 1). The series in the elliptic nome for the form factors contributing to C(0, 1) 
were first conjectured in sec. 5.2 of [4]. At that time it was not clear how to sum up such 
expressions^. Here we have used a completely different method to obtain a closed form expression 
for C(0, 1) in terms of elliptic functions. 

Following [24] we will use two variables s and t; 

4 j sm\i{2E /keT) , T > T^, 

t = s , s = < 1 (4) 

[smh{2E /ksTy, T <Tc. 

Let us introduce a special notation for diagonal correlation functions 

C^it) = C{N,N), (5) 

where ± denotes T regimes, respectively. We note in the symmetric case the dualof C^(i) 



^The idea of a A-extension was made even earlier in McCoy, Tracy and Wu [19] in connection with solutions 
to Painleve III. 

^ After this paper appeared, the technique for summing up such series was given by McCoy et al. in [23]. 
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Our main result is the following Toda-type relation for diagonal correlations 

t-^ logC^(t) + - logC^(t) + = ^^—^ (c^(,))2 - (6) 

valid for both temperature regimes at N = 1,2,.... Supplemented by initial conditions for 

iV = 0, 1, relation (6) gives a very effective way to calculate C^{t) for iV = 2, 3, 4, It took us 

under 10 mins to calculate Cj^{t) for TV = 2, . . . , 20 on a laptop computer using a Mathematica 
program. 

What's more important is that eqn (6) remains valid for the A-extension of diagonal corre- 
lation functions. In fact, as shown in section 5, a form factor expansion for A?^ = is sufficient 
to calculate the generalized diagonal correlations for = 1, 2, . . .. 

The paper is organized as follows. In section 2 we define the A - extended pair correlation 
functions, give integral formulae for the form factors and recall the result of Jimbo and Miwa 
on diagonal correlators. In section 3 we list some basic facts on the Hamiltonian formulation 
of Painleve VI theory and prove a Toda-type relation for tau-functions related by a special 
shift transformation. In section 4 we establish a connection with the Jimbo-Miwa result and 
prove relation (6). In section 5 we calculate generalized diagonal correlators at arbitrary A for 
iV = 0, 1, 2. In section 6 we present a closed form expression for the A - extension of the first 
nontrivial off-diagonal correlation function M = 0, = 1. Together with the results of Section 4, 
these provide all initial conditions required for the quadratic difference equation satisfied by the 
A-extended correlation functions, and so allows, in principle, the rapid calculation of any required 
correlation function. The major advance is in the calculation of A-extended correlation functions. 
In section 7 we discuss the results obtained, and mention possible further developments. 

2 Form factor expansions of pair correlation functions 

We start with the A - expansion of pair correlation functions in multi-particle components 

oo 

C+(M,A^;A) = s-^(l-t)V^^A2"+^C'(2"+^)(M,A^), T > n, (7) 

n=0 

and 

oo 

C-(M,A^;A) = (l-t)V4(i + ^A2"C'(2")(M,iV)), T < T^, (8) 

n=l 

where for n > (see, for example, formula (4.2) in [4], given here for the symmetric case) 



C(")(M, N) 



TT TT 

1 f duji f duJn 

'n\J ~2tt"'J '2tF 

— TT —IT 



n 



L sinh7i 



^ 2 

n 

l<i<j<n 



n 

cos(iV^a;i), (9) 



i=l 



Sinh7i = (y2 _ 1)1/2^ = - (yf - l)^^ ^ yi = S + S ^-COSUi, (10) 

We shall call C^(M, A^;A) generalized correlation functions. The pair correlations for the 
symmetric Ising model are obtained by setting A = 1 

C^{M,N) = C^{M,N;l). (12) 
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Using a different approacfi Jimbo and Miwa [18] in 1980 introduced tlie function 



t(i-l)^logC+(iV,iV)-i T>n, 
t{t-l)j^logC-{N,N)-^, T<T,. 



(13) 



and showed that it satisfies the "sigma" form of Painleve VI 

(Pa\'^ .da f , ^da 1\ f da \ ,^9/, ^da 



In 2004 [25] Forrester and Witte gave a new derivation of this result working from the Toephtz 
determinant form rather than the Predholm form. 

As mentioned in the introduction, Jimbo and Miwa also considered an isomonodromic A- 
extension of C^{N,N) which also satisfies (14). Using extensive computer calculations Boukraa 
et al. [20] observed that the Jimbo-Miwa A-extension is the same as the form factor expansions 
(7-8). It provides an opportunity to investigate properties of generalized diagonal correlation 
functions ^^{NjN-jX) by studying solutions of (14). 

In the rest of the paper wc will assume A to be arbitrary unless its value is explicitly speci- 
fied. We shall also use the same notation introduced in (5) for generalized diagonal correlation 
functions (7-8) 

C^{t) = C^{N,N;\) (15) 

and simply omit the dependence on A. 

The solutions of (14) can be specified by their expansion near t — t- 0. In [20] Boukraa et al. 
analyzed series expansions of C^(t) at t ^ and demonstrated that 

C-^t) = (1 - t)'/' + >^'^^§ff§^t^^'i^ + Oit)) (16) 

and 



^2 



cut) = A(i - t)y^fS^, + >^'^^0§§f^t^''^'^\^ + om (17) 

where 

/J;.';, = <'"^^2F,(i,Ar + i,Af + l;t), (18) 

and {a)N = r(a + N)/T{a) is the usual Pochhammer symbol. We shall use (16-17) to fix a 
normalization of recurrence relations in section 4. 



3 Some facts from Painleve VI theory 

In this section we shall briefly review the properties of the Hamiltonian form of the Painleve VI 
(Pvi) equation [26,27], given by eqn (22). Our main reference is the work of Okamoto [28]. 

For each solution of the Pyi equation, one can construct a Hamiltonian function H{t) = 
Hit; b) which depends on four parameters 

h = {huh2MM)- (19) 
A tau-function r(t) = r(t; b) related to the Hamiltonian H{t; b) is defined by 

i?(i;b) = ^logr(t;b). (20) 
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Let us introduce an auxiliary Hamiltonian h{t) = h{t;h), 

h{t) = t{t - l)H{t) + 62(61, 63, 64) t - ^62(61, 62, 63, 64) (21) 

where ei{xi, . . . is the i-th. elementary symmetric function in n variables and a set of Xj's 
can be a subset of 6j's as in (21). 

The function h{t) solves the equation Ey7[b]^ 

h'{t) t{l - t)h"{t) ^ + h'{t)[2h{t) - {2t - l)h'{t)] + 61626364 ^ = (^'(*) + ^k) ■ (22) 

fe=i 

The group G of Backlund transformations of Py/ is isomorphic to the affine Weyl group of 
type F4: Wa{F4). It contains the following transformations of parameters (not all of them are 
independent) 

lui : 61 -H- 62, tU2 : 62 -H- 63, i(;3 : 63 -H- 64, 104 : 63 ->■ —63, 64 —64, (23) 
: kq Ki, : kq Kqo, x^ : ^ 6 (24) 

where 

/«o = 6i + 62, Ki = 61 - 62, too = 63 - 64, 6* = 63 + 64 + 1 (25) 
and the parallel transformation 

/3 : b = (61, 62, 63, 64) ^ b+ = (61, 62, 63 + 1, 64). (26) 

For each transformation s from the group Wa(-^4) such that 

s : b H^- bs (27) 

one can construct another auxiliary Hamiltonian h{t\hs) which satisfies Ey/fbg]. The function 
h{t;hs) is a rational function of t and 

h{t-h), JtHt;h), ^h{t;h)Y (28) 

If we consider a sequence of transformations related to /™, m = 0, ±1, ±2, . . . and construct an 
associated sequence hjn{t) and related tau-functions Tm{t), then it is well known [28] that such 
tau-functions satisfy Toda-type recurrence relations. 
Let us consider a different transformation Ig G G: 

Is = X1X2I3X2X1, (29) 

where the action of the Xj's and ^3 on 6i's is defined by (23-26). 
After simple calculations one obtains 

If^ ■ (61, 62, 63, 64) ^ (61 ± ^, 62 T ^,63 ± ^, ^4 ± ^). (30) 

To calculate the expressions for related auxiliary Hamiltonians we need to use the results 
from Okamoto's paper [28], sections 2,3. Here we shall only present the answer 



^So named by Okamoto [28] in his studies of PVI. It is in fact a generalization of the "sigma" form of PVI 
(14). 
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h±{t) = h{t, b±) = h{t) - - ± + 



1 ^ Si 2t{t - lfh"it) ± [25i/i(t) + 253(1 -t)- S1S2 



2{2{t - l)h'{t) - 2h{t) + S2) 



where 



Si = 6,(61,-62,63,64), b± = lt\h) = (61 ± ^,62 T ^,63 ±l,h± I). 



Now we can follow the Okamoto derivation of Toda relations for tau-functions. Using 
and (31) we obtain 



h+{t) + h^{t) - 2h(t) = t{t - 1) 



H+{t) + H^it)-2H{t) + 



3t 



+ 



Substituting 



2 2{t - l)h'{t) - 2h{t) + ^2 



H^{t) = I log[r±(t)], H{t) = I log[T(t)] 



and integrating over t we get 



,T+(t)r_(t) 



B{t-l)^^^^Ll_yLl = it-l)h'{t)-hit) + '' 



S2 
2 



or after using the expression for h{t) in terms of r(t), 



d'^ , /X d ^ , , „ T+(t)T-(t) 

where B is the integration constant. 

We have just shown that if we introduce a sequence of parameters 



On= 1 61 + —,62 - y,63 + + y 



and associated sequences of Hamiltonians and tau-functions 

hN{t) = hit,hN), HN{t) = H{t,hN), TN{t) = T{t,hN), 

then tau-functions rjv(t) satisfy Toda- type recurrence relations 

d^ ... d TN+lit)TN-l{t) 
'dfi + dt = ^) ' 



where the normalization constant can depend on the parameters 61,62,63,64 and A''. 



4 Recurrence relations for C^{t) 

Now we use the results from the previous section to obtain recurrence relations for generalized 
correlation functions C]^{t) valid for arbitrary A. 

First we need to establish a connection between equations (14) and (22). Let us consider 
the equation Ey/fbAr] (22) with parameters (37) and choose 

6i = -l/2, 62 = -1/2, 63 = 0, 64 = 0. (40) 

If we denote a corresponding solution of (22) as /liv(t), then after a substitution 

1 

hN{t) = aN{t)-—t + -, (41) 

equation (22) coincides with (14). 

Now combining formulae (13, 20-21) and (41) we can find a connection between tau-functions 
r7v(t) from the previous section and C^(t), viz. 



Substituting (42) into (39) we finally arrive at recurrence relations for generalized diagonal 
correlation functions C^(t), 

t^^,logC^(t) + ^^logC^(t)+^^_^^, - ^^^^^^^^ . (43) 

The normalization constant Bn in (39) can be calculated by substituting expansions (16-17) at 
t — > into (43) . It appears that 

Bn = N^-^ (44) 

for both the T > Tc and T <Tc regimes. 

Let us comment on these Toda-type relations for C^{t) which are valid for arbitrary values 
of A. They allow us to calculate diagonal correlations for any = 2, 3, 4, . . . starting with Co{t) 
and Ci{t). We don't need to use correlations at the dual temperature s^^ as in [18]. In fact, 
when A is arbitrary, it is sufficient to know C^{t) only, as C^{t) can then be calculated from 
(43) as we shall see in the next section. 

To our knowledge, equations (43) have not previously appeared in the literature. They 
provide a very efficient method for the calculation of diagonal form-factors and sums in formulae 
(7-8) in closed form. 



5 The Picard solution and calculation of Cg (^) and C^{t) 

In this section we shall calculate initial conditions for the recursion (43). The functions Cf{t) 
are simply related to tau-functions ro(i) by (42). A related function ho{t), given by eqn (41) 
solves the equation Ey/fb] (22) with parameters bi given by eqn (40). It is well known that a 
general solution of (22) depending on two arbitrary parameters (specified by initial conditions) 
can be found in this case. 



8 



Consider another set of parameters 6j 



61 = 0, 62 = 0, &3 = -^, h = -l. (45) 

This choice corresponds to the Picard solution [29] of the Painleve VI equation which can be 
found in closed form (see, for example, [30]) and expressed in terms of the elliptic Weierstrass 
function depending on two arbitrary parameters. 

One can represent solutions of Painleve equations in terms of logarithmic derivatives of ratios 
of functions with no movable singularities. As one of many examples of such a representation, 
we mention a representation of solutions of the Pyj as logarithmic derivatives of tau-functions 
(see section 4.3 of [28]). 

A location of zeros of tau-functions for the Picard solution (and related the Hitchin solution) 
has been investigated by Kitaev and Korotkin [31] and recently by Brezhnev [32]. It follows 
that the tau-function for the Picard case should be proportional to the Jacobi theta-function 
depending on parameters specified by initial conditions. A detailed calculation of tau-functions 
for the Picard class solutions has been given recently by the first author in [33]. Here we will 
use these results. 

First, let us denote an auxiliary Hamiltonian and related tau-function for the Picard case (45) 
as hp{t) and Tp(t) respectively. Since the equation (22) is symmetric in the fej's we immediately 
conclude that 

ho{t) = hp{t). (46) 
Using (21) and integrating over t we produce a relation between tau-functions 

ro(i) = (l-t)VVp(t). (47) 

We introduce an elliptic parameterization and define the elliptic modulus k as 

k = s^, t = k^. (48) 

Following [34] we will define theta functions Oi(x\T) with quasi-periods tt and ttt. The elliptic 
modulus k in (48) and its complement k' are 

mv?' "^3(o|r)2' ^ ^^^^ 

It will be more convenient to use the parameter t as the second argument of the elliptic functions 
^ = ^ = K(t) = |^i(0|T), K'(i) = i^(l-t), (50) 

where K(i) and K'(i) are the complete elliptic integrals of the first kind of parameters i and 
\ — t respectively. 

Now let us introduce Jacobi elliptic functions 

where 

u = — (52) 

TT 
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and define the fundamental elliptic integral of the second kind [34] as 

S{u,t) = / dn^ {x,t)dx. (53) 
Jo 

It satisfies 

S{K{t),t) = E{t), (54) 

where E(f) is the complete ehiptic integral of the second kind. 

It is convenient to use slightly renormalized complete elliptic integrals 

A- = HM, E^ml. (55) 

TT TT 

One of the results of [33] can be written as 

rHt) = c.(x,»)//''t-'/*^l|ii^. (56) 

where x and y are two arbitrary parameters and co(x, y) is a normalization constant independent 
oft. 

First let us consider the high temperature regime T > Tc- Combining (42, 47, 56) we obtain: 

C.,).«(...)<i^//"?i^ (57) 

for some values of x and y. Recall that C(^(t) depends on the parameter A. It remains to find 
a relationship between the parameters x and y on the RHS and the parameter A on the LHS 
of (57). We can expand both parts of (57) in series in x, y and A and compare the expansions 
term by term. 

Integrating the well known formula for a logarithmic derivative of 94,{x\t) we obtain 

xK 

e4(x|r) = e4(0|r)exp{-yi^^+ j S{y,t)dyY (58) 



Using (51) and (58) we can expand the rhs of (57) in a series in x and yr with polynomial 
coefficients depending on t, K and E. 

For C(^(t) we need to calculate the first few form factors in the sum (7). To that end, we can 
use the integral representation (9) or a different integral representation [20] for the functions 

/g"/i)=(7P"+i)(7V,7V)/., f'^^^^ = C^^^\N,N). (59) 

It was observed in [20] that the diagonal form factors /^J^jy are polynomials in t, K and E 
for small values of n, and the first four, corresponding to n = 1 ... 4 are given in [20]. Now 
we can expand the integral representation for the diagonal form-factors in series in t up to 
sufficiently high order and compare it with the above polynomial Ansatz. In that way one can 
produce closed form expressions for the form-factors. The first few exact form factors C'''"')(0, 0), 
n = 1, 2, . . . , 7 are given in Appendix A. We note that in [4], eqn (5.7), Orrick et al. gave a 
conjectured expression for C^"^ (0, 0) giving the results in terms of q series. Evaluation of their 
expression for n = 1,2, ... ,7 gives a series expansion that agrees with the expansion of the 
expressions given in Appendix A, thus providing mutual conformation. 
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Now matching series on both sides of (57) up to order we achieve perfect agreement if we 
choose 



Q 7 

tLi tJU Jb ^ / Q \ 



co(x,y) = l, A = x-- + — -^ + 0(x^), y = 0. (60) 



Prom this expansion we can guess that 



and we obtain 



A = sin(a;) (61) 



CjW = C-(0,0;A)=(i^M=2g^, (62) 

a result which has not been published before. A nontrivial self-check of the relation (61) is to 
set A = 1 which corresponds to the symmetric Ising model. 
We obtain 

C (0,0,1)- ^4(0|r) - tV4 eMr)~ ^ ^ 

as expected. 

Another comment is in order. Based on connections to Painleve VI theory we have proved 
that the generalized diagonal correlations can be represented in terms of theta- functions. How- 
ever, the relations between the parameters x, y and A are still conjectural. One of the possible 
ways to prove such relations is to find a partial differential equation for generalized correlations 
which includes derivatives with respect to both A and t. 

We can handle the low-temperature regime T < Tc in a similar way. The only difference is 
that now we have different relations between parameters, 

co{x,y) = —iexp{ix), A = sin(a;), y = 7r/2. (64) 

The answer reads 



. _ , . _ , ,,i/4 ^4(arcsin(A)|T) _ g4(arcsin(A)|r) 

Co(t)-C (0,0,A)-(l-t) ^-^^ - ^-^^ 



(65) 



Again for A = 1 

Co~(i) = C-(0,0;l) = Mi^ = l. (66) 

The formula (65) first appeared in Boukraa et al. [20] (see formula (92), Section 6). The 
authors of [20] remark that they have calculated A-extended diagonal correlations for N = 
0, 1, 2, . . ., and have given a selection of the algebraic curves which arise from special choices of 
A.4 

We emphasize that in our derivation of (62, 65) we relied on the exact calculation of tau- 
functions for the Picard solution [33]. This provides an alternative derivation of equation (65), 
previously announced in [20]. 

As mentioned above we have conjectured relations between parameters (60-61) and (64) 
based on series expansions in A and x and the Ising case A = 1. However, formulae (63) and 
(65) produce analytic expressions for the integrals (9) at M = = for any positive integer n. 
We checked numerically up to n = 10 that we obtain correct values for the form factors in (9). 



* After this paper appeared, results for Af = 0, 1 diagonal correlations were given in [23], which follow from 
conjectured results in [4]. 
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Now we turn to the calculation of C^{t). Let us start with the auxiliary Hamiltonian ho{t) 
corresponding to the Picard solution. We take the tau- function (56) with y = (y = tt/2) and 
arbitrary a; for T > Tc {T < Tc) and use eqns. (20-21) with parameters (40) to find /io(t). Now 
we can find auxiliary Hamiltonians h±{t) corresponding to AT = ±1 from (31-32). It is easy to 
see that they coincide 

h+{t) = h-{t). (67) 
It follows from (20-21) and (37) with N = ±1 that 

Ti{t) =a{x)tT-i{t) (68) 

and from (42) 

Cf(t) = a±{x)C%{t) (69) 

for some integration constants a±{x). 

Substituting into (43) at A'" = we obtain 

We can calculate first and second derivatives of Cq (t) . For that we need the following differen- 
tiation formulae 

#«=^--. (71) 

dt 2t{l-t) 2t dt 2t ^ ^ 

d , , sn('U, t)cn^('U, cn(u,t)dn(u,t) r , ^, 
^sn(„, ,) = - '^(il) + W-1) - 1) + ^(«- "1 ■ 

|.n(„,<)= sn^(M)cn(-. _ «) _ ^ ^,„_ ,,3, 

cZ , , , sn^(ti, t)dn(u, t) sn(tt, t)cn(n, t) r , ^, ,t 

c/„ _ £!(M)_ ^,„2/ , sn(^,^)cn(u,t)dn(^,^) 

^ " 2(i^ ~ 2 2(t^:T) ^^^^ 

and 

d_l f e4{x\T) \-] _ sn^xK,t) {S{xK,t)-xEf 

V^4(o|T)yJ 4(1 -t) 4t(t-i) ■ ^ ^ 

The last formula can be proved using (58) and a partial differential equation for theta-functions. 
After straightforward simplifications we arrive at 

t^\ogC^{t) + ^logCo-(t) = -4^^^ (cn(z,t)dn(z,t) +sn(z,t)X)2 (77) 

and 

4.ogC+(*) + |logC^W = -j^^ (78) 

where 

z = xK ='^^^^^, X = ^[\ogei{x\T)]'=£{xK,t)-xE. (79) 
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Substituting back into (70) we have 

1 e4{x\T) 



and 



VM^^3(0|r) 



C+(t) 



(cn(z, t)dn{z, t) + sn{z, t) X) 



i{x\t) X 



(80) 



(81) 



VM^^3(0|r) tV2- 

To determine the factors a±{x) we use (61), expand (80-81) in A and compare with form factor 
expansions (7-8) at M = A'^ = 1. These form-factors are given in [20] up to order nine in A and 
we will not reproduce them here. All these expansions are consistent with a choice 



^/cL\Jxj = -s/ a-{x) = cos(x). 



(82) 



Again wc can do an independent check by expanding near the point .t = 7r/2 or A = 1 
which corresponds to the symmetric Ising model. In this case we obtain the results for the 
nearest-neighbor diagonal correlations C(l,l) 



C+(l,l;l) =i-i/2 ^^t-l)K + E] 



and 



C-{1,1;1)=E, 



(83) 
(84) 



which agree with formulae (4.9a) and (4.9b), Section VIII. 4 of [11]. 

Now we can substitute the initial values of C^(t), C^{t) into recurrence relation (43) and 
calculate C]^{t) for any A?^ = 2, 3, ... for arbitrary A. All we need are the differentiation formulae 
given above. 

We obtain for = 2 and T <Tc 



C2{t) 



ei{x\T) 



cos2(a;)6'3(0|T) 



+ -C S X + X"^ ( 2S^ 



t + l 

3t 



X^ 
3t 



(85) 



and for T > Tr 



c^t) = - 



U(x\r) 



3t cos2(x)6i3(0|r 
+ S[X^ + X^{6tS^ -5{t + l)] +tC^ 



2CX [2X^ + 2tS'^ - t-l] + 



where X is defined in (79), and 



C = cn{z,t) dn{z,t), S = sn{z,t). 



(86) 



(87) 



The answer for C^{t) is given in Appendix B. 

Let us point out that the factor C^_^(t) always divides the Ihs of (43). We don't have a 
rigorous proof of this cancelation. Exactly the same phenomenon occurs for many polynomial 
tau-functions related to Painleve equations. Usually the proof of this fact relies on determinant 
representations of such tau-functions. It would be extremely interesting to find a determinant 
representation of C^(i) for arbitrary A. 

As mentioned in the Introduction, the equations (43) for A = 1 together with initial con- 
ditions (63, 66, 83-84) and differentiation formulae (71) give a very efficient procedure for the 
calculation of diagonal correlation functions {N, N) for the symmetric Ising model. 
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6 A calculation of C±(0, 1; A) 

In this section we perform a resummation of the A-extension of the simplest row nearest- 
neighbor correlation C='=(0, 1;A). It is not known whether ofF-diagonal correlations and their 
A-extcnsions satisfy nonlinear differential equations. Therefore, a knowledge of explicit expres- 
sions for C^{M,N;X) at small values of M, N can be useful. 

First we shall start with the expansion of integrals (9) at s — )• 0. We have 

sinh7j = — scoswj + 0(s^)), Xi = ^(1 + scoswj + 0(s^)), 

hij = sin^ {~2~^^ + cosLOj) + 0(s^)). 

Then it is easy to see that C^'^\0, 1) has a kinematic zero of order tt,^ as s ^ 0. 

Expanding C^") (0, 1) in s for n = 1, 2, 3, 4, 5, 6 one can see that C'(")(0,1) has a zero in s of 
order n(n + 1) and guess the first two nonzero terms of the expansion 

C'(")(0,l) = ^^^(l + ^ + O(54)) at s^O, n = 1,2,3,... (89) 

We expect these integrals to be polynomials in K, E and [20]. The case n = 1 is easy 

Ca)(0,l) = -i + l(l + .2)i^, (90) 

and the answer for n = 2 was given in [20] 

^^(')(0, l) = l-\{l + s')K - i(.2 - 3)(1 + .2)^2 _ (91) 

Using series expansions in s up to a sufficiently high order we were able to obtain similar 
expressions for n = 3, 4. They are given in Appendix B. The structure of the n = 1, 2, 3, 4 cases 
suggests the following Ansatz for general n : 

(7^(0, 1) = af\EKy + (1 + .2) ^ ^ 6g)£;^K-V,_2,-(.) (92) 

1=0 i=l 3=0 

where [x] stands for the integer part of x, o!f^ and 6^"^ are rational coefficients, and Pn{s) is the 
polynomial in of degree n — 1 with integer coefficients: 

Pn{s) = Co,„ + Ci,„s2 + . . . + ^_2^^g2(n-2) ^ ^2(n-l)_ (93) 

It is quite remarkable that the condition for the sum in (92) to have a zero of order n(n + 1) 
in s fixes all coefficients in (92) together with polynomials Pn{s) up to a common normalization 
factor. This factor is determined by the coefficient of the leading term in (89). 

Applying this procedure we easily calculated C^^\0, 1) for n = 5, 6, . . . , 15. The coefficients 
Qj^^ and b^^'j together with polynomials Pn{s) up to n = 6 are listed in Appendix B. 

Of course, this is just experimental mathematics based on matching series expansions. It is 
desirable to run some self-consistency tests. We evaluated the integrals C^"\0, 1) forn = 1, . . . , 6 
numerically at s = 0.95 and they matched (92) with 1 part in 10^* accuracy which corresponds 
to approximately 360 terms in the series expansions in s. So the Ansatz (92) is almost certainly 
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correct. Note that eqn (5.8) in Orrick et al. [4] gives a conjectured expression for ^^(O,!) 
in terms of Jacobi theta functions and q series. Expanding that equation in a (5'-scrics and 
converting the series expansions of (Bl) to (B7) to g-series provides a valuable consistency 
check. We confirm that the series agree to the first 100 terms. 

Now having explicit expressions for the form-factors we can try to perform a resummation 
of the next-neighbor correlations C(0, 1; A). 

As we did for the diagonal correlations, it is natural to look at the the logarithms of 
C^(0, 1; A) and expand them in series in x, where 

X = arcsin(A). (94) 

Then we obtain 

KF 

log C±(0, 1; x) = ^— x2 + /±(x, K, s^) (95) 

where functions f^{x,K^s'^) as series in x do not depend on the elliptic integral of the second 
kind E. The absence of E in the functions f^ix^K^s^) is a very strong argument in favour of 
our Ansatz (92). 

Using definitions (7-8) and (58) we can write 

^^+(0,1;.:) = \j^^F^i^^K,s% C-{0,l;x) = {x, K, s'). (96) 

Now let us rewrite scries expansions for F'^i^x, K, ,s^) in variables x and y = xK and split them 
into even and odd parts in y. It appears that the resulting two functions factorize and have 
similar structure 

F±(x, K, s') = a±{x) [gevenixK, s^) - p^{x) 9oddixK, s^)] . (97) 

where a±{x) and P±{x) depend on the variable x only and geveniv, s^), Qoddiy^ ■s^) do not depend 
on X. A similar separation happens for generalized diagonal correlations at A?^ = 0, 1 and is 
another strong argument in favour of our calculations. 

Using Ansatz (92) we calculated form factors (9) for n = 1, . . . , 15 and obtained series 
expansions for functions a±{x) and P±{x) up to x^^. Then it is easy to guess that 

Finally we have to find the functions geven{y, s'^) and godd{y,s^) where the variable y = xK. 
Their expansions in y have polynomial coefficients in of growing complexity and it is not easy 
to reconstruct them. 

However, based on the diagonal case we assumed that these functions are double-periodic in 
y and substituted their expansions into the differential equation 

4 

{dyg{y,s)f = Y,A{s){g{y,s)y, (99) 

i=0 

where the coefficients Ai{s) are independent of the variable y. 

This leads to immediate success, and we can easily find the polynomial coefficients Ai{s) and 
solve the differential equation (99) for both cases. 
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Collecting all the above formulae we get the final answer 



where 



C-(0,l;x) 

c+{o,i;x) 



eA{x\q) cos(f) 
^3(01?) cos(x) 

s 63(01^) cos(a;) 



cn(2^,s4)dn(^,s4) 



gev{x, s) - tan^^l^ godd{x, s) 
gev{x,s) - cot(^|^ godd{x,s) 



godd{x, s) 



(l + s2)gn(^|^,s4) 

l + s2sn(^,s4)2 



(100) 



(101) 



l-s2sn(^,s^)^ 

and the variable x is related to A by (61) or (94). 

The final check is to expand (100) near x = Tr/2 which corresponds to the symmetric Ising 
model and compare the answer with known nearest-neighbor row pair correlation functions. We 
obtain 



and 



C+(0,1) 



c-(o,i) 



2s 



(102) 



(103) 



which match (4.5a) and (4.5b) of Section VIII.4, [11]. 

Let us briefly comment on the formulae (100). To our knowledge closed expressions for 
generalized diagonal correlations C^(0, 1) have not appeared in the literature before. 

As discussed in in the introduction in 1980 McCoy, Perk and Wu obtained quadratic differ- 
ence relations for correlation functions. In our notation the A-extension of these equations is 
given in [4]. These difference equations allow one to calculate generalized correlation functions 
provided that wc know the answer for diagonal correlations and the first off-diagonal pair correla- 
tions C(0, 1). Therefore, formulae (100) provide the initial conditions for a recursive calculation 
of all generalized pair correlation functions C(M, N; A) for arbitrary M, N = 0,1,2, .. . 

While undoubtedly correct, a rigorous proof is still lacking. The missing ingredient is, for 
example, a Toda type relation similar to that obtained above for the diagonal correlations, 
plus verification of the connection between the parameters A and x. This is currently under 
investigation. 

After the appearance of the first draft of this paper, it was pointed out to us by W.P. 
Orrick and J-M. Maillard that the result we give in eqn (100) is, in principle, derivable from eqn 
(5.8) in the paper [4] by Orrick et al. by summing over the correlation functions and making 
appropriate substitutions. We have not confirmed this result. We remark however that if that 
is indeed the case, further exact expressions for C(1,1;A), C(2, 0;A), and C(2, 1;A) should be 
similarly derivable from eqns (5.9)-(5.11) of [4] and consistent with Perk's difference equations. 



7 Conclusions and outlook 

In this paper we have found and proved a set of Toda-type difference equations (6) satisfied 
by diagonal pair correlation functions. It is interesting that the A- extension of these pair 
correlations satisfy the same equations. Wc think that there should be a direct way to prove 
relations (6) from the Toeplitz determinant representation of the diagonal correlations. That 
would then enable us to understand the determinantal representation for their A-extension. 

When that is done, we could try to generalize relations (6) for the case of row correlations 
where the Toeplitz determinant representation is known for A = 1. However, it could well be 
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difficult to generafize it to arbitrary values of A. In fact, formulae (100) provide the first examples 
of closed form expressions for the A- extended nearest-neighbor row correlation functions. We 
failed to find a nonlinear differential equation satisfied by (100) valid for any value of A. Of 
course, that doesn't mean that such an equation doesn't exist. 

Let us emphasize that based on connections to Painleve VI theory we have proved the 
results in equations (57, 65), up to the non-rigorous identification of the connection between the 
parameter A and the parameters x and y. 

We also note that in 2005 Forrester and Witte [35] gave a nonlinear system of recurrences in 
their corollary (6.3), related to the discrete Painleve V equation for the same ratio of correlations 
that appears on the rhs of (43). 

Finally, as we mentioned above, the derivation of the generalized diagonal correlations and 
relations (100) is not rigorous, based as it is on series expansions. However it is undoubtedly 
correct, satisfying numerous checks, and should provide initial conditions to calculate all gen- 
eralized pair correlations C{M,N;X) for any M and N from the quadratic difference equations 
of Perk. This should provide all necessary tools for a significantly more efficient generation of 
the series expansion for the A-extended susceptibility. An ambitious, but still distant hope is 
to sum the correlation function series and produce the long sought closed-form expression for 
the susceptibility. We are planning to report on the evaluation of C{M,N;X) in a forthcoming 
publication. 
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Appendix A. The diagonal case 

In this Appendix we present the first seven form factors (43) at M = = used for the 
calculation of initial conditions for the recursion relations (6). We also give generalized diagonal 
correlation functions for N = 3. 
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Form factors C'(")(0, 0) forn = 1, . . . , 7 : 
sK 

K ~ T ~ 2"' 

C'(^)(0,0)_K tK^ E K'^E KE'^ 

K ~ 6" ~8 VT ~ 6" ^ ^ 8 



(7(6) (0,0) 



K 



I5" ^ T2" 



+^1-45 + 



(- 


2t\ 

3"y ^ 48" 




KE 


K^E^ 




12 


48 


6 



22t 8t2 

~ Ts" Ts" 

+ —tt;- + 



16 



+ 



16 



(7(7) (0,0) 



sK 



5 
112 



+ 



240 



48 



+ {f-U + Q)K^{ — 



K^E^ 
48~ 
KE\ 



it - 2)K' — 



37 



5KE K^E^ 
+ 



72 



48 



+ 



144 240 J 5040 



{t - 2){f - lot + 10)K^ 



(A.2) 



In [36] it is pointed out by Bostan et al. that the minimal order ODE for the n-particle 
contributions to the susceptibihty %(") are achieved by appropriate Unear combinations of %("') 
with terms x^^~'^\ X^'^^^K etc. In particular, these combinations arc: x(^) — x^^^/(i, X^"^^ "X'^^V^j 
-^(5) _-^(3) /2+x^^^ /120, x(^^ — 2x(^V3+X*'^'' /45. It is remarkable that the analogous combinations 
of form factors C'(")(0, 0) for n = 1, . . . , 7 give expressions that are homogeneous in products of 
elliptic integrals of degree n. ^ That is to say, combinations such as (7(6) (0,0) — 2(7(^)(0, 0)/3 + 
(7(2) (0, 0)/45 contain only terms of the form K^, K^E, K^E^ and K^E^. Prom (A.2) it is easy 
to conjecture that the appropriate linear combination that will give the minimum order ODE 

for x^^^ is X^^^ - 5x^^V3! + 13x(^V5! - X^^V^! 

Generalized diagonal correlations C^{t) are obtained by solving (6) for N = 2. They are: 



and 



135*5/2 cos^{x)93{0\t) 

04{x\t) 



135*2 cos^{x)e3{0\T) 



tC SRl + XR+ 

C + S 



(A.3) 

(A.4) 



This was pointed out to us by B G Nickel. 
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R+{t) = 168X6 + 84X^ [6tS^ -7{t + l)] +U{l-S^){l-tS^){2tS'^-9{t + l)) 
+108X2 [2tS'^{tS'^ -3{t + l)) + 3 + At + 3t'^] , (A.5) 



Rt{t) = -X^ + 6X^[6tS^-t-l]+ 9X'^ [42 1 S^{t - t - 1) - + 28t - l] 

+ 2X^ [l05t'^S'^{2tS^ - 5{t + 1)) + 3*5^(99 + 196 1 + 99 1^) -{t + l){2 + 103t + 2*^)] 

+ 3t{l- S^){l-tS^) [7*52(3*5^ - ll(t + l)) + 24*2 + 21t + 24] ^ ^^ g^ 

R^{t) = 9X^ + 42X^ [2*5^ - t - 1] + 3X^ [14*5^(3*52 - t - 1) - 13*^ + 28* - 13] 

+ 6X2 ^3 ^2^4,^2 * 52 + * + 1) - 3 * 52(3 + 52* + 3*2) + (* + l)(2 + 13* + 2*2)] 

+ *2(l-52)(l-*52) [*54 + 3 52(* + l)-135] , (A.7) 

R- (t) = X^ + eX^ [6 * 52 - 5(* + 1)] + 9X^ [14 * 52 (* 52 - * - 1) + *2 + 28 * + 1] 

+ 2X2 t'^s'^{2S^t - * - 1) - 3 * 52(13 + 84* + 13*2) + (* + 1)(20 + 61 * + 20*2)] 

+ 3*(l-52)(l-*52) [*52(3*52 + 5(* + l))-4*2-197*-4] . (A.8) 



Appendix B. The off-diagonal case 

In this Appendix we first give form factors C'^"^(0, 1) for n = 1,2,3,4 obtained by series 
expansions of (9) to series in s up to 200 terms. 

CW(0,1) = -1 + 1(1 + .2)^, (B.l) 

(B.2) 

C'('H0,1) = ^ - 1(1 + .2)i^-l(.2_3)(l + , 2)^2 (B_3) 

(B.4) 

^(^^(0,1) = + 1^(1 + .2)K + 1(.2-3)(1 + . 2)^2- (B.5) 

(B.6) 

C^'hO,!) = ^s-fQi^ + s')K-^{s'-3){l + s^)K^+ (B.7) 

+^is^ + 6s2 - 11)(1 + s2)x3 + 1 (56 _ 21^4 ^ 3^2 ^ 25)(1 + s2)K^ - 
96 oo4 

-—EK + -E'^K'^ + -EK'^{1 + s2) + —EK^{s'^ - 3)(1 + ^2) 
48 8 8 16 

We note that, in analogy to the hnear combinations of form factors ^^"^(0, 0) for n = 1, . . . , 7 
that give expressions that arc homogeneous in products of elliptic integrals of degree n reported 
in Appendix A, rather more complicated combinations of form factors C^") (0, 1) produce similar 
results for the form factors C^'"\0, 1). In particular, we note that 

C(3)(o, 1) + lc(2)(o, 1) - ^C'(i)(0, 1) - ^ = 0{EK',K') 
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n 


(n) 
% 


(n) 


(") 


(n) 


6 


231 


3319 


25 


1 


1024 


11520 


192 


48 


5 


63 


23 


1 




256 


96 


16 




4 


35 


17 


1 




128 


48 


8 




3 


5 

16 


1 
4 






2 


3 
8 


1 
2 






1 


1 
2 









Table 1: Coefficients for n = 1, . . . , 6. 



n 


An) 
"1,0 


An) 
"2,0 


"3,0 


^(n) 
"3,1 


^(n) 
"4,0 


An) 
"4,1 


'^5,0 


^(n) 
"5,1 


"5,2 


An) 
"6,0 


"6,1 


^(n) 
"6,2 


6 


1289 


3319 


3 


9 


25 


25 


1 


1 


1 


1 


1 


1 


7680 


46080 


256 


64 


9216 


384 


7680 


192 


32 


46080 


768 


64 


5 


263 
1280 


23 

384 


7 
384 


7 

32 


1 
768 


1 

32 


1 

3840 


1 

96 


1 

16 








4 


19 
96 


17 

192 


1 

96 


1 

8 


1 
384 


1 

16 














3 


13 
48 


1 

16 


1 

48 


1 
4 


















2 


1 

4 


1 
8 






















1 


1 

2 

























Table 2: Coefficients 6)7 forn = 1, . . . , 6 



and 

C(^)(0,1) + \C^'\0,1) - f^C('H0,l) - ^CW(0,1) + ^ = 0{EK^E'K\K% 

Here we also list coefficients a,["^ , b^-"j and polynomials Pn{s) for n = 1, ... ,6 which occur in 
eqn (92) for the form factors C^''\0, 1). 



Pi{s) = 1, 

P2{s) = -3 + s\ 

P3{s) = -11 + 6s^ + s\ 

P4{s) = 25 + 3s^ - 21s'' + s^, 

P5{s) = 201 - 180s^ - 66s'' + 60s^ + s^ 

pe{s) = -299 - 123s2 + 466s'' + 106s^ - 183s^ + s^°. (B.8) 
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